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form  a  vapor  jet  which  could  then  serve  as  the  conducting  channel  in  a 
giant  power  pulse  discharge. 

A  model  based  on  an  energy  balance  Is  presented  which  predicts  the 
mass  removed  and  the  pressure,  density  and  temperature  of  the  material 
as  it  propagates  from  the  target  surface  as  functions  of  the  pulse  time 
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and  the  absorbed  power  density.  Combinations  of  laser  pulse  time  and 
power  density  are  considered  for  which  the  surface  temperature  remains 
low  enough  In  order  to  avoid  cut-off  by  a  plasma  plume  In  front  of  the 
target  surface. 

A  model  of  adiabatic  jet  expansion  Is  presented  to  predict  the  radial 
expansion  of  the  vapor  plume  as  It  propagates  In  the  axial  direction. 
Temperature,  pressure  and  density  are  predicted  after  radial  expansion 
using  the  model  of  adiabatic  jet  expansion. 

These  models  predict,  theoretically,  that  the  desired  masses  and 
limits  on  vapor  plume  expansion  can  be  achieved  using  relatively  low 
power  densities  for  the  proposed  experiment,  — 
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ABSTRACT 


This  study  is  to  assess  the  feasibility  of  using  a  laser's 
radiation  to  remove  masses  on  the  order  of  milligrams  from 
a  target  surface  and  to  form  a  vapor  jet  which  could  then 
serve  as  the  conducting  channel  in  a  giant  power  pulse  dis¬ 
charge  . 

A  model  based  on  an  energy  balance  is  presented  which 
predicts  the  mass  removed  and  the  pressure,  density  and 
temperature  of  the  material  as  it  propagates  from  the  target 
surface  as  functions  of  the  pulse  time  and  the  absorbed  power 
density.  Combinations  of  laser  pulse  time  and  power  density 
are  considered  for  which  the  surface  temperature  remains  low 
enough  in  order  to  avoid  cut-off  by  a  plasma  plume  in  front 
of  the  target  surface. 

A  model  of  adiabatic  jet  expansion  is  presented  to  predict 
the  radial  expansion  of  the  vapor  plume  as  it  propagates  in 
the  axial  direction.  Temperature,  pressure  and  density  are 
predicted  after  radial  expansion  using  the  model  of  adiabatic 
jet  expansion. 

These  models  predict,  theoretically,  that  the  desired 
masses  and  limits  on  vapor  plume  expansion  can  be  achieved 
using  relatively  low  power  densities  for  the  proposed  experi¬ 


ment. 
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I .  INTRODUCTION 


Requirements  for  material  and  equipment  tests  at  very 
high  x-ray  radiation  fluxes  with  particular  spectral  charac¬ 
teristics  have  led  to  the  development  of  giant  pulse  power 
devices,  Johnson  [1]. 

The  technique  of  using  "exploding  wires"  or  wire  arrays 
dates  back  to  1773  [2] .  In  this  technique  the  energy  from 
a  large  storage  capacitor  is  discharged  first  through  a  pulse 
forming  network  then  through  thin  metallic  wires  which  pinch 
due  to  the  development  of  large  magnetic  fields  and  eventually 
develop  into  a  high  temperature  plasma  which  becomes  a  source 
of  high  intensity  x-rays.  In  recent  time  the  wire  array  has 
been  replaced  by  a  gas-puff  from  a  supersonic  nozzle.  After 
preionization  with  microwave  radiation  the  discharge  through 
the  gas  jet  leads  to  a  high  temperature  pinch  discharge. 

The  x-ray  intensity  and  spectral  distribution  depends 
upon  the  high  voltage  pulse  form,  plasma  temperature  and 
dynamics,  mass  of  the  wire  array  or  gas  in  the  gas-puff, 
z-number  of  the  material  and  the  overall  geometry  of  the 
experiment.  Elaborate  computer  codes,  Colombant  [3],  have 
been  written  to  model  the  x-ray  emission  as  a  function  of  these 
parameters.  In  order  to  obtain  the  desired  spectra  and  inten¬ 
sities,  these  computer  codes  predict  masses  on  the  order  of 
milligrams  and  high  z-number  materials  are  required  in  the 
discharge  channel. 


Wire  arrays  are  not  available  in  small  enough  diameter 
sizes  for  high  z -number  materials.  As  for  the  gas-puff 
device,  materials  again  are  severely  limited  to  gaseous  ele¬ 
ments  such  as  Ar,  Kr,  and  Xe.  Because  of  the  shortcomings 
of  wire  arrays  to  satisfy  the  mass  and  z-number  requirements, 
alternate  methods  of  establishing  high  energy  discharge 
channels  are  of  interest. 

The  use  of  high  power  laser  pulses  directed  at  a  target 
material  has  been  proposed  by  Schwirzke,  described  by  Johnson 
[1] .  In  such  a  scheme,  the  target  material  in  the  form  of  a 
solid  disc  can  be  mounted  on  the  anode  (or  cathode)  of  a  high 
voltage  pulsed  power  machine.  A  laser  beam  is  directed  toward 
the  target  and  a  certain  mass  of  the  target  is  vaporized, 
figure  1.  This  vapor  plume  will  expand  through  the  evacuated 
chamber  toward  the  cathode.  Electrical  energy  would  then 
couple  to  the  vapor  plume  producing  a  pinched  discharge 
channel  with  an  intense  x-ray  emission.  Preionization  of  the 
vapor  could  be  done  by  microwaves  as  in  the  gas-puff  experi¬ 
ment  or  conceivably  by  the  laser  itself.  The  advantages  of 
this  concept  over  wire  arrays  and  gas-puffs  are  great  in 
principle  as  targets  could  conceivably  be  made  of  almost  any 
material. 

By  controlling  the  energy,  pulse  width  and  focal  area  of 
the  laser  along  with  the  choice  of  material,  the  x-ray  emission 
spectrum  could  be  varied.  Combinations  of  different  materials 
could  also  be  used. 


There  are  two  distinct  mechanisms  by  which  a  laser  cam 
remove  material  from  a  target.  First,  the  laser  energy  is 
absorbed  by  the  target  to  increase  the  temperature  to  the 
vaporization  temperature.  At  high  power  densities  this  takes 
a  very  short  time,  on  the  order  of  a  few  nanoseconds.  Since 
the  power  density  is  large,  a  high  percentage  of  ionization 
in  the  vapor  cloud  or  plume  will  be  reached  after  a  very  short 
time.  This  high  ionization  level  causes  absorption  of  the 
laser  energy  in  the  plume  thus  shielding  the  target  surface 
from  the  continuing  incoming  radiation.  However,  further 
vaporization  of  the  surface  via  electronic  heat  conduction 
continues.  This  mechanism  has  been  investigated  by  Brooks 
[4].  He  used  a  Q-switched  neodymium  glass  laser  with  a  5 
Joule  output  of  pulse  duration  25  nsec  directed  at  an  aluminum 
target  with  a  focal  diameter  of  0.15  cm.  This  results  in  a 
power  density  of  about  1010  Watts/cm2.  Measurements  showed 
mass  removal  to  be  about  0.7  mg.  While  the  mass  yield  is  of 
the  desired  order,  an  important  feature  of  this  mechanism  makes 
it  less  attractive  for  the  purpose  here.  The  main  plasma, 
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carrying  the  bulk  of  the  ionized  mass,  propagated  at  1.1x10 
cm/sec.  If  anode-cathode  separation  were  assumed  to  be  4  cm, 
the  main  plasma  would  reach  the  anode  in  364  nsec.  The  target 
temperature  is  so  high  that  vaporization  continues  for  almost 
10  iisec  after  the  laser  pulse  is  stopped.  As  a  result,  the 
plasma  particles  in  the  plume  will  expand,  thus  filling  the 
entire  cavity  before  sufficient  mass  has  been  removed  from  the 
target.  As  the  plasma  now  covers  a  large  volume,  energy 


coupling  from  the  high  voltage  source  will  become  difficult. 
Contamination  of  and  arcing  in  the  gap  of  the  coaxial  delay 
line  forming  the  cathode  cone  will  reduce  efficiency. 

An  alternate  approach  may  be  used  to  overcome  this  prob¬ 
lem.  If  the  laser  power  density  is  lowered  the  ionization 
level  can  be  reduced  to  a  point  where  no  significant  absorp¬ 
tion  of  the  laser's  radiation  takes  place  in  the  vapor  plume. 
This  plume  is  now  a  true  vapor  plume  and  not  a  plasma.  The 
laser  energy  is  not  shielded  from  the  target  surface.  As  the 
surface  temperature  is  now  relatively  low,  the  plume  propagates 
much  slower.  The  timing  of  the  application  of  the  high  voltage 
can  be  made  to  correspond  to  the  travel  time  of  the  vapor 
plume.  The  question  is  whether  significant  masses,  on  the 
order  of  milligrams,  can  be  removed  from  the  target.  This  is 
the  approach  that  will  be  studied  in  this  paper. 

There  are  many  parameters  controlling  the  mass  removal 
process.  The  energy  delivered  by  the  laser  to  the  material 
surface  divides  into  the  reflected  portion  and  the  absorbed 
portion.  Classically  the  absorption  is  dependent  on  the  laser 
radiation  wavelength  and  the  conductivity  of  the  target  mater¬ 
ial.  In  section  II  ,  the  relevant  absorption  processes  are 
described  for  various  materials  of  interest  and  a  table  of 
absorption  is  compiled  which  is  used  as  the  basis  for  the 
absorption  calculated. 

The  absorbed  energy  at  the  surface  raises  the  temperature 
of  the  material  leading  eventually  to  melting  and  vaporization 
beginning  at  the  surface.  The  initially  steep  temperature 
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gradient  near  the  surface  leads  to  rapid  heat  conduction 
into  the  material.  This  process  is  governed  by  the  heat 
conduction  equation.  Because  of  the  three  phases  present 
this  is  a  very  complex  problem.  There  are  many  models  of 
this  process  using  various  approximations  in  the  literature, 
often  employing  elaborate  computer  programs.  In  order  to 
stay  with  the  purpose  of  this  study,  to  assess  the  general 
feasibility  of  the  proposed  mechanism,  a  very  simple  model  of 
the  heat  conduction  and  evaporation  has  been  adopted  which 
will  be  discussed  in  Section  III.  Its  relation  to  other  models 
in  the  literature  are  also  discussed. 

As  the  material  evaporates  from  the  surface  it  expands 
both  axially  and  radially.  Section  IV  deals  with  the 
radial  expansion  which  is  described  as  an  adiabatic  jet  expan¬ 
sion.  Density,  temperature  and  pressure  changes  in  the  expan¬ 
sion  process  are  also  described. 
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II.  THE  ABSORPTION  PROCESS 


Materials  can  be  divided  into  three  basic  classes, 
conductors,  semi-conductors  and  dielectrics.  Dielectric 
materials  are  not  of  interest  here  because  they  do  not  con¬ 
duct  well.  In  order  to  vaporize  these  materials,  high  power 
densities  must  be  used  which  would  result  in  plasma  production 
at  the  target  surface.  This  is  a  condition  we  are  trying 
to  avoid  as  these  high  temperature  plasmas  lead  to  surface 
shielding. 

Semi-conducting  materials  in  general  are  also  poor  conduc¬ 
tors.  However,  at  certain  wavelengths,  there  are  resonances 
where  conduction  is  quite  high.  In  these  regions  their  thermal 
properties  resemble  those  of  metals.  Therefore,  our  considera¬ 
tions  will  be  confined  to  the  metallic  absorption  process. 

The  interaction  of  materials  with  electromagnetic  waves 
is  formulated  on  the  basis  of  Maxwell's  equations.  This 
development  is  well  known  and  documented  [5,  6,  7].  Derivation 
of  optical  quantities  is  presented  in  Appendix  A. 

Plots  can  be  made  of  absorbtivity ,  A  vs.  X,  the  wavelength, 
for  several  materials  using  the  following  expression: 

A  *  - ^ - 7  (eq.  2-1) 

(n+ir  +  k 

where  n  and  k  are  the  real  and  imaginary  parts  of  the  complex 
index  of  refraction.  The  n  and  k  are  both  functions  of 


frequency,  w,  the  conductivity,  a,  and  the  parameter  m*/N 
where  m*  is  the  effective  mass  and  N  is  the  free  electron 
density  of  the  metal.  In  general  it  is  very  difficult  to 
obtain  an  appropriate  value  for  m*.  Much  controversy  has 
arisen  on  how  to  measure  effective  mass,  Seitz  [8]  lists  some 
approximate  values  which  are  shown  in  Table  1.  Harrison  and 
Webb  [9]  review  analytical  and  experimental  work  which  illus¬ 
trates  the  difficulty  in  obtaining  appropriate  values.  This 
problem  of  choosing  a  value  for  m*  need  not  present  an  obsta¬ 
cle,  however.  Wieting  and  Schriempf  [7]  show  proof  that  the 
sensitivity  of  reflectivity  to  the  parameter  0/N,  where 
0  *  m*/mQ,  is  small.  Here  mQ  is  the  free  electron  mass. 
Therefore,  the  selection  of  values  for  m*  in  [8]  are  justified 
for  our  purpose. 

Figures  2-9  are  plots  of  n,  k,  and  A  as  functions  of 
decreasing  wavelength  using  the  data  from  Table  I.  It  is 
interesting  that  the  values  for  absorbtivity  are  very  low  for 
the  region  where  w  is  smaller  than  the  plasma  frequency.  The 
absorbtivity  then  increases  very  sharply  at  the  plasma  fre¬ 
quency  for  the  electrons  in  the  metal  and  quickly  reaches  100 
percent.  The  degree  of  absorption  remains  relatively  constant 
over  the  infrared  region  and  is  small,  generally  1-2  percent. 

The  question  now  arises  as  to  whether  the  calculated 
values  for  absorbtivity  are  reasonable.  Bennett,  Silver  and 
Ashley  [10]  have  shown  surprisingly  close  correlations  between 
experiment  and  these  classically  calculated  absorbtivities  in 
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ultra  high  vacuums  of  10”^  torr  for  aluminum.  However,  at 

lower  values  of  vacuum  the  theory  breaks  down  and  the  absorb- 

tivity  is  greater.  This  is  significant  because  in  this  experi- 
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ment  the  base  pressure  in  the  vacuum  chamber  is  ~10  torr. 
Bonch-Bruevich  [13]  also  shows  experimental  results  illus¬ 
trating  the  change  in  reflectivity  as  a  function  of  pulse  time 
in  copper  and  steel,  figure  10.  It  must  be  pointed  out,  how¬ 
ever,  that  these  results  pertain  to  conditions  different  from 
those  of  interest  here.  This  particular  experiment  was  carried 

out  using  a  laser  with  a  power  density  on  the  order  of  10^ 

2 

Watts/cm  .  At  this  high  power  level  there  may  be  major  plasma 
production  and  the  marked  decrease  in  reflectivity  at  the  tar¬ 
get  surface  is  probably  due  more  to  the  plume  absorption  rather 
than  the  temperature  dependent  changing  of  optical  properties 
of  the  target  material.  The  first  plateau  at  2  ysec  is  probably 
due  to  the  onset  of  vaporization.  One  would  not  expect  a  much 
lower  power  level  laser  to  produce  the  degree  of  change  in 
absorbtivity  that  Bonch-Bruevich  observed. 

In  calculating  the  absorbtivity,  dc  conductivity  was  used. 
Conductivity  is  inversely  proportional  to  the  resistivity. 

The  resistivity  of  metals  increases  nearly  linearly  with 
temperature  in  the  region  from  room  temperature  to  melting. 
Figure  11  is  a  plot  of  absorbtivity  as  a  function  of  temperature 
for  four  metals,  aluminum,  silver,  cadmium  and  zinc,  all  at 
X  *  10.6  pm.  The  underlying  assumption  in  this  plot  is  that 
the  liquid  and  solid  phase  conductivities  are  the  same.  The 
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FIGURE  11.  ABSORBTIVITY  vs.  TEMPERATURE  AT  10 


plot  shows  increases  in  absorbtivity  on  the  order  of  1-3 
percent. 

Thomas  and  Musal  [14]  have  proposed  using  an  effective 
value  for  absorbtivity,  A^ff,  using  the  following  expression: 

Aeff  “  (Av-VIln(W  (eq.  2-2) 

where  Av  is  the  absorbtivity  at  vaporization  temperature  and 
Aq  is  at  room  temperature.  Harrach  [15]  points  out  that  these 
effective  values  of  absorbtivity  may  be  5-10  times  those  at 
room  temperature  and  cites  the  Bonch-Bruevich  report,  [13] , 
as  experimental  evidence.  Again,  this  could  be  misleading 
at  low  power  densities. 

Surface  preparation  plays  a  large  role  in  the  absorbtivity 
of  any  practical  metal  surface  and  is  an  empirical  matter. 
Schriempf  [16]  states  some  values  for  Al,  Au,  Cu,  Ag,  Ti, 

Mg  and  304-steel.  Sandblasting  the  surface  appears  to  be  the 
best  technique  resulting  in  up  to  ten  times  the  abosrbtivi- 
ties  of  ideal  surfaces.  He  does  not  state  the  pressure  under 
which  his  values  for  absorbtivity  were  found. 

In  this  paper  four  materials  have  been  chosen  to  illus¬ 
trate  various  relationships.  Aluminum  is  used  because  of  the 
vast  amount  of  previous  efforts  in  the  laser  weapons  program. 
Silver,  cadmium  and  zinc  were  chosen  because  of  the  wide 
variation  in  their  thermal  properties. 


III.  THERMAL  RES PONSES 


A.  THE  MODEL 

After  the  laser  energy  strikes  the  target  surface,  the 
absorbed  energy  is  converted  into  thermal  energy  in  the 
material.  The  thermal  energy  so  deposited  in  the  surface 
is  distributed  in  various  ways  in  the  material.  The  immedi¬ 
ate  result  is  a  sharp  rise  of  the  surface  temperature.  This 
causes  a  steep  temperature  gradient  in  the  material  which  leads 
to  heat  transport  from  the  surface  and  into  the  material  via 
heat  conduction.  Eventually  the  surface  region  reaches  melting 
temperature  and  finally  vaporization  temperature  and  the  sur¬ 
face  melts  and  vaporizes  when  enough  energy  has  been  absorbed 
to  overcome  the  heat  of  melting  and  vaporization.  The  exact 
development  of  this  3  phase  process  and  the  temperature  profile 
into  the  medium  as  a  function  of  time  is  governed  by  the  heat 
conduction  equation  together  with  appropriate  boundary  condi¬ 
tions  . 

The  heat  conduction  equation  in  three  dimensions  is: 

pg  Cg  =  div(K  grad  T)  +  A  (eq.  3-1) 

where  p  is  a  solid  density,  C  is  the  heat  capacity,  K  is 

9  9 

the  thermal  conductivity  and  A  is  the  heat  produced  per  unit 
volume  per  unit  time  by  the  laser  source. 


A 


(1  -R)l 


4irk  -4nkz/X 


3F 
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(eq.  3-2) 


where  (1  -  R)  is  the  absorbtivity,  I  is  the  laser  power  density 
2 

in  Watts/cm  and  F  is  the  absorbed  power  density.  The  term 
4irk/A  has  the  units  of  cm  *  and  is  defined  as  the  absorption 
coefficient  where  k  is  the  imaginary  part  of  the  complex  index 
of  refraction.  The  value  S  *  A/4irk  is  the  skin  depth.  The 
skin  depth  is  that  distance  where  the  intensity  of  the  laser 
is  reduced  to  1/e  of  its  value  at  the  target  surface.  In 
metals,  this  value  is  on  the  order  of  10  ®  cm  for  wavelenghts 
in  the  infrared. 

There  are  different  boundary  conditions  for  each  phase 
of  the  target  material.  In  general  even  the  source  term.  A, 
is  very  complicated  as  absorbtivity  varies  with  temperature 
and  phase.  The  laser  pulse  is  directed  at  the  target  at  an 
angle  and  the  pulse  shape  is  Gaussian  in  nature.  The  target 
thickness  effects  the  temperature  profile  in  the  target  mater¬ 
ial.  Even  the  material  constants,  C_,  p.  and  K  change  with 

s  s 

the  wide  temperature  ranges  and  phase  changes  involved  between 
room  temperature  and  vaporization  temperature. 

In  its  full  generality  the  problem  is  quite  intractable 
and  for  the  purpose  here,  significant  simplifications  must  be 
made  by  adopting  a  number  of  assumptions  and  approximations 
which  will  now  be  presented  and  discussed. 

We  will  begin  the  discussion  by  considering  the  development 
of  the  temperature  profile  for  a  single  solid  phase  material 
until  vaporization  temperature  is  reached  at  the  surface,  that 
is,  we  neglect  the  complication  of  the  liquid  phase.  Specifically 
we  make  these  assumptions: 
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1.  It  will  be  assumed  that  the  laser  pulse  is  constant 
with  time  and  the  absorbtivity  is  also  constant,  possibly 
an  effective  value  calculated  by  the  method  of  Thomas  and 
Musal  [14]  previously  mentioned. 

2.  The  entire  problem  will  be  done  in  one  dimension. 

This  one  dimensional  approach  is  valid  if  the  laser's  beam  is 
perfectly  uniform  over  a  large  focal  area. 

3.  The  materials  thermal  parameters,  p  ,  C  ,  K,  etc., 

s  s 

are  assumed  to  be  temperature  independent. 

4.  Starting  temperature  for  the  process  will  be  assumed 
to  be  at  T  =  0°  so  that  all  temperatures  will  be  given  in  °C. 

5.  The  model  considered  is  a  laser  beam  directed  at  a 
semi-infinite  solid  target  material,  figure  12.  Under  these 
conditions  an  exact  solution  of  the  heat  conduction  equation 
can  be  obtained,  Carslow  and  Jaeger  [17].  It  is  also  assumed 
that  the  laser  wavelength  is  in  the  infrared  region  where  the 
skin  depth  is  very  small,  on  the  order  of  10  6  cm. 

Under  these  conditions  the  heat  conduction  equation 
reduces  to 


32T  1  3T  _  -(1  -  R)I  ct  e 

77  "  <  3t  *  K. 


(eq.  3-3) 


where  <  ■  K/p  C  an<*  a  the  absorption  coefficient  4irk/X. 
s  s 

The  initial  and  boundary  conditions  are  T  *  0  at  z  =  o  and 


K 


3T 

II 


(1  -  R)I  -  F  at  z  *  0. 
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FIGURE  12.  MODEL  SHOWING  LASER  RADIATION  DIRECTED  AT 
A  SEMI- INFINITE  SOLID  AND  THE  CREATION  OF 
THE  VAPOR  FRONT  IN  DISTANCE  Az 
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The  solution  to  the  problem  is: 


T  (z ,  t ) 


2F  _ 

— /ict/ir 


e  z  /^<t  -  ^  erfc  z/2/ict} 


(eq.  3-4) 


where 


erfc  =  1  -  erf 


and 

o  x  •  2 

erf  -  —  /  e"x  dx' 

✓?  o 

This  model  allows  us  to  estimate  the  time  when  the  surface 
reaches  vaporization  temperature. 


2F  _ 

T(0,tv)  =  Tv  »  /7t7¥ 


(eq.  3-5) 


or 


v2m  2 
K  T 

ir  v 

*  7T1 


(eq.  3-6) 


Braginskii,  et  al.,  [18],  shows  experimental  support  for 
this  equation.  Harrach  [19]  uses  the  so-called  "heat  balance 
integral"  method  to  solve  the  problem  of  moving  boundaries 
and  obtains  a  similar  result  for  tv: 


tv  -  (3  [1  -2e’1]/KHKTv/Fo}2  (eq.  3-7) 

The  numerical  ratio  of  the  two  solutions  for  t  is  0.99. 
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The  difficulty  with  this  solution  to  the  heat  conduction 
equation  is  that  it  does  not  yet  involve  phase  changes.  In 
order  to  remove  mass  from  the  target  surface,  first  the 
vaporization  temperature  must  be  reached,  i.e.,  enough  energy 
must  be  deposited  into  the  target  surface  to  raise  the  tem¬ 
perature  from  Tg  »  0  to  Tg  *  Tv,  or  greater  (Tg  =  surface 
temperature).  In  addition,  enough  energy  must  also  be  de¬ 
posited  to  overcome  the  latent  heats  of  fusion  and  vaproiza- 
tion.  We  now  need  a  simple  model  for  this  second  stage  in 
the  process. 

Figure  12  illustrates  the  model  of  a  semi-infinite  solid 
material.  As  vaporization  is  reached,  there  will  be  an 
instantaneous  removal  of  material  from  the  target  surface. 
This  results  in  a  moving  boundary  condition.  One  would  now 
have  to  solve  a  heat  conduction  problem  with  these  moving 
boundaries,  where  the  solid  interface  (again  neglecting  the 
liquid  phase)  moves  at  a  recession  velocity  U  . 

If  this  velocity  can  be  calculated,  then  one  could  deter¬ 
mine  the  mass  removal  rate.  The  solution  of  this  heat  con¬ 
duction  problem  with  its  moving  boundary  in  one  dimension  has 
been  given  by  Harrach  (15]  in  analytical  form  using  again  the 
"heat  balance  integral"  method.  An  estimate  of  this  velocity 
can  be  obtained  from  a  simple  heat  balance  expression.  If 
Az  is  the  distance  that  the  interface  moves  in  time  At  then, 


,  £t 

o  Az 


L  P_  +  c_  T  p_ 

s  S  V  3 


(eq.  3-8) 


28 


where  L  is  the  total  latent  heat  from  room  temperature  to 
vaporization  temperature,  L  -  1^  +  Lv.  Cg  is  some  overall 
heat  capacity  and  p_  is  an  overall  density  (C  and  p„  may 
be  assumed  to  be  those  values  for  the  solid  to  simplify  the 
problem) .  Then : 


VL  +  csV 


(eq.  3-9) 


(eq.  3-10) 


This  velocity  could  be  thought  of  as  a  steady  state  velocity. 
This  expression  for  the  steady  state  recession  velocity  of 
the  evaporating  surface  is  identical  to  the  one  Harrach  ob¬ 
tains  from  his  model.  However,  in  Harrach' s  model,  at  times 
close  to  the  vaporization  time,  the  surface  recession  velocity 
is  lower.  As  no  closed  form  expression  exists  for  times  near 
the  vaporization  time,  he  shows  a  plot  obtained  numerically 
to  illustrate  the  buildup  of  the  recession  velocity.  This 
solution  is  reproduced  in  Figure  13.  The  result  of  this 
analysis  is  that  the  value  for  the  steady  state  surface 
recession  velocity  obtained  from  the  general  heat  balance  must 
be  reduced  in  the  region  where  total  laser  pulse  time  is  on 
the  order  of  surface  vaporization  time. 

The  considerations  discussed  above  do  not  take  into  account 
yet  that  the  material  immediately  on  the  target  surface  may 
be  at  a  temperature  higher  than  the  vaporization  temperature. 
Anisimov  [20]  suggests  that  the  surface  is  indeed  at  a  higher 


Normalized  recession  velocity 


V°s 

steady  state 


Normalited  time  t(t/ty)-l] 


FIGURE  13.  NORMALIZED  TIME  vs.  NORMALIZED  SURFACE 
RECESSION  VELOCITY.  NUMERICAL  SOLUTION 
FROM  REF.  16. 
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temperature  caused  by  an  elevated  boiling  point  created  by 
the  dynamic  and  high  pressure  situation  set  up  at  the  target 
surface.  We  will  now  expand  our  simple  model  somewhat  so 
that  it  incorporates  the  idea  that  the  surface  may  be  at  a 
higher  temperature  than  the  vaporization  temperature  and 
that  the  vapor,  when  it  does  leave  the  surface,  leaves  at 
a  velocity  of  U  *  /k_T  /m.  This  would  be  the  z-direction 

V  o  S 

2  2  2  2 

component  of  the  rms  thermal  velocity  where  vx  +  v  +  vz  *  <v> 
and  kfi  is  the  Boltzman  constant.  We  make  the  following 
assumptions  for  this  new  model: 

B.  THE  THEORY  OF  EVAPORATION 
1.  Assumptions 

We  will  now  assume  that  with  the  laser  energy  directed 
at  a  target  an  amount  of  material  of  depth  Az  is  raised  above 
the  vaporization  temperature  and  then  is  removed,  Figure  12. 

The  materials  thermal  properties,  as  before,  are 
assumed  to  be  constant  with  temperature.  The  liquid  and  solid 
phase  densities  and  heat  capacities  are  considered  to  be  approxi¬ 
mately  equal  and  the  gas  phase  heat  capacity,  Cv,  is  equal 
to  3/5  C  .  These  assumptions  are  based  on  data  obtained  by 
experiment,  Kelly  [21],  CRC  Handbook  [11],  and  JANAF  Thermo¬ 
chemical  Tables  [22] .  In  addition,  it  is  assumed  that  the 
material  in  the  small  increment  az,  to  be  removed,  has  heat 
capacity  and  solid  density,  pg,  as  it  has  not  had  time 
to  expand  significantly  while  being  heated  to  a  temperature 


31 


2.  Temperature  Profile  in  the  Material 

Under  the  above  conditions  and  the  one  dimension 
model  of  semi-infinite  target  thickness,  one  would  expect 
the  temperature  profile  to  be  as  before  (eq.  2-4): 


T(z,t)  =  ^{£t7^  e"z2/4Kt  -  |  erfc [z/2/xt]  } 


This  expression  is  cumbersome  to  work  with  as  it  contains 

the  error  function.  However,  it  can  be  approximated  fairly 

closely.  Let  D,  the  diffusion  length,  be  the  distance  where 

the  temperature  decays  to  1/e  of  its  surface  value,  T  . 

s 

Solving  the  above  equation,  the  distance  D  is  found  to  be 

D  ■  0.97/Kt  or  simple  D  z.  /let.  It  is  readily  apparent  that 

the  temperature  profile  is  of  a  decaying  exponential  form. 

If  the  diffusion  length  is  used  as  the  spatial  decay  constant, 

-z/D 

the  equation  T  *  Tge  can  be  compared  to  the  actual  solution. 
The  comparison.  Figure  14,  illustrates  the  close  approxima¬ 
tion.  In  fact,  the  difference  is  only  a  few  percent  and  even 
less  if  one  is  interested  in  a  small  distance,  Az  from  the 
surface.  Figure  15. 

3.  Energy  Balance 

An  energy  balance  equation  can  be  written  to  describe 
the  heat  conduction  problem  using  the  temperature  profile 
T  »  t  e-z/D: 


00  Z 

FQt  ■  psLAz  + psCs  /  T8e”z/Ddz  + Psc-TvAz  + pscv  /  (Tae“z/D  -T  )dz 

Az  0 


TEMPERATURE  (°C) 


2400 


FIGURE  14.  COMPARISON  OF  ACTUAL  SOLUTION  AND  APPROXIMATE 
SOLUTION  FOR  T(z)  IN  THE  SEMI-INFINITE  SOLID 
MODEL 
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D 

DISTANCE  AND  TARGET 
ARBITRARY  UNITS 


TEMPERATURE  PROFILE  IN  A  SEMI-INFINITE  SOLID 
SHOWING  THE  SMALL  VAPOR  LAYER  Az 


In  this  equation,  the  terms  pgLAz  and  PgcgTvAz  are  the  energy 
required  to  raise  the  small  increment  Az  to  vaporization 
temperature  and  overcome  the  latent  heats.  The  term 

ps°s  0/42<v'2/d  -  Vd2 

is  the  energy  required  to  raise  Az  to  a  temperature  Tg  >  Tv* 
The  second  term, 

PsC,  /"v‘*/D  dz 

Az 

is  the  energy  lost  due  to  heat  transport  into  the  material. 
This  balance  holds  in  particular  for  t  *  tp,  the  laser  pulse 
length.  Using  the  substitutions 


T  « 

V 

Ta  e“Az/D 

(eq. 

3-12) 

D  = 

(eq. 

3-13) 

Az  = 

D  ill  (Tg/Ty) 

(eq. 

3-14) 

Cv  - 

3/5  Cs 

(eq. 

3-15) 

the  above  energy  balance  equation  reduces  to: 


The  parameter  FQ/t^  may  now  be  used  to  calculate  the  surface 
temperature  from  this  energy  equation.  Implicit  in  this 
derivation  is  that  the  pulse  duration,  t  *  t  is  comparatively 
small  while  Fq  is  large  so  that  the  vaporized  material  has  no 
time  to  escape  during  the  duration  of  the  pulse  and  so  all 
the  energy  arriving  must  be  absorbed  in  the  same  material 
leading  to  a  temperature  rise  above  the  vaporization  tempera¬ 
ture. 

Alternatively,  keeping  FQ/t^  constant  one  could  pro¬ 
duce  a  laser  pulse  where  t  *  t  is  very  long  and  Fq  small. 

In  that  case  the  vaporized  material  escapes  from  the  surface 
before  the  pulse  terminates.  The  laser  energy  heats  the  next 
surface  layer  to  the  vaporization  temperature  and  as  it  is 
reached  the  surface  layer  escapes  and  so  on,  so  that  the  sur¬ 
face  temperature  stays  essentially  at  the  vaporization  tem¬ 
perature.  The  equation  Az  *  D  ln(Tg/Tv)  obviously  breaks  down 
in  this  large  tp  case  as  it  predicts  Az  =  0.  The  equation 

for  F  t  in  this  situation  is 
o  p 


F  t 
o  P 


PgLAZ  + 


P  C 
^s  s 


Az 


► 


0,c,Tv“ 


(eq.  3-17) 


The  expression  can  now  be  used  to  calculate  Az//t ~  as  a 

tr 

function  of  FQ/t^,  namely! 


Fo^  -  ".Vv^ 

VL  +  csV 


(eq.  3-18) 
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4.  The  Mass  Removed 

The  mass  removed  during  the  total  time,  tp,  can  be 
expressed  as  M=p  Az  where  M  is  the  mass  expressed  in  grams 
per  laser  focal  area.  In  the  case  of  short  pulses  we  can 
substitute  Az  «  J>ln  (T  /T  )  . 


M  *  pg  /<tp  ln(Tg/Tv) 


(eq.  3~19) 


we  now  express  M  as  a  function  of  the  parameter  FQ/t^. 


MFq  *  FQ/t^  pg  ln(Tg/Tv)  (eq.  3-20) 


It  must  be  kept  in  mind  that  Tg  is  also  a  function  of  FQ/t^. 

In  the  case  of  long  pulses,  when  the  surface  tempera¬ 
ture  stays  close  to  vaporization  temperature  we  use  the 
expression  for  Az//t^  given  above. 


MF. 


PsAzFo 


FOv%1  (L+CgTv) 


]  (eq.  3-21) 


5.  Axial  Expansion  of  the  Vapor  Jet 

For  the  purpose  of  the  intended  application  it  is  of 
importance  to  estimate  the  distance  which  the  vapor  front 
travels  during  the  duration  of  the  laser  pulse.  If  that  dis¬ 
tance  is  small,  on  the  order  of  Az  or  less,  then  no  significant 
expansion  has  taken  place  during  the  energy  deposition  and  we 
can  use  the  mass  removal  relation  of  the  high  surface  temperature 
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model  where  T  >  T  .  In  this  case  there  will  be  no  further 
s  v 

experimental  constraint.  One  would  fire  the  main  discharge 
of  the  high  voltage  source  when  the  vapor  jet  has  expanded 
axially  across  the  gap  which  would  be  long  after  the  end  of 
the  laser  pulse. 

In  the  case  of  long  pulses  it  becomes  important  how 

far  the  vapor  surface  has  traveled  when  the  laser  turns  off 

and  all  mass  has  been  vaporized.  That  distance  traveled 

must  be  commensurate  with  the  experimental  geometry  to  give  a 

feasible  operating  condition. 

If  t  is  the  total  pulse  time  and  t  is  the  time  to 
p  v 

vaporization  then  (t  - 1  )  *  tf  is  the  excess  time  remaining 
after  vaporization  in  which  the  mass  located  in  z  can  be 
removed  from  the  target  and  travel  a  distance,  1,  from  the 
target  surface.  The  vaporized  particles  will  travel  at  a 
velocity  approximately  equal  to  the  rms  value  of  velocity 
in  a  one  dimensional  Maxwellian  distribution. 

Uv  *  /kBTv/m  (eq.  3-22) 

Schriempf  [16]  has  obtained  the  same  result  for  vapor  velocity 
in  slightly  different  form  using  simple  heat  balance  equations 
and  conservation  of  mass: 


/rt7z 


(eq.  3-22') 


where  Z  is  the  atomic  weight  of  the  target  and  R  is  the  gas 
constant,  R  *  8.314  x  10^  erg  °K  mol  The  time  to  vapori¬ 
zation  may  be  calculated  when  the  expression  (eq.  3-18)  is 
applied  to  the  time  tv  when  the  surface  reaches  vaporization 
temperature,  T  =  T  . 

S  v 


t 


v 


(eq.  3-23) 


using  again  K  =  <p  C  .  This  value  is  fairly  consistent  with 

s  s 

the  previously  stated  values  (eq.  3-6)  obtained  by  Carslow 
and  Jaeger  [17]  and  Harrach  [15]  using  completely  different 
methods,  and  with  the  experimental  results  of  Braginskii  [18]. 
An  expression  for  the  excess  time,  tf,  can  now  be  written: 


t 


f 


„2  _  2 

K  T 

fcp - T* 

r  <  F_ 


(eq.  3-24) 


The  distance  which  the  vapor  surface  travels  in  the  axial 
direction  during  the  pulse  time  is  then 


i  =  Uy(tp  -  tv)  (eq.  3-25) 

We  will  consider  a  "short"  pulse  when  l  £  Az.  For 
the  criterion  for  determining  whether  we  are  in  the  long  or 
short  pulse  region,  we  can  safely  use  the  expression  for  the 
long  pulse  which  does  not  involve  the  elevated  surface 
temperature. 
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p_/tT  -  p_c  t  /< 
41  *  ^ -» 


Combining  the  above  expression  for  8,  and  Az  we  have  « 
for  short  pulse: 


K2T  2  F  /t~  -  p  C  T 


which  can  be  rewritten  to  give  the  condition. 


Fo  ± 


F  Jt~[Fo-/tp,"  PsCsTv^'  ] 
°  PC  P.Tt  +  CaTv> 


If  this  condition  is  not  satisfied  we  must  calculate 
value  of  l  where 


1  -  wv 


K2  T  2 
Uv(tp  “  < - 7~ 


or 


U 


•  ■ 

Fo 


6.  Density  and  Pressure  in  the  Plume 

Once  expressions  are  found  for  Az  in  the  two 


(Eq.  3-26) 


condition 


(eq.  3-27) 


(eq.  3-28) 


the 


(eq.  3-29) 


(eq.  3-30) 


cases, 


long  and  short  pulse,  then  we  can  calculate  the  density  and 
pressure  of  the  vapor  plume  of  the  surface.  Conservation  of 


mass  tells  us  that 


P  U 
v  v 


psUs 


(eq.  3-31) 


where  pv  is  the  vapor  density  as  Az  expands  into  the  cavity. 
In  the  case  of  short  pulses  Az  *  D  ln(Tg/Tv),  where  D  */ict  ', 
Ua  =  A^/(tn-tv)  and  Uv  *  /m.  Substituting,  one  can  find 


p  v' 


the  expression  for  py. 


ps^  ln(VV 

Vtp  -  tv) 


v2  _  2 
K  Tv 


(eq.  3-32) 


Pv  =  Pg  FQ/t;[^  ln(Ts/Tv)] 
F0  Uv  [(F0/t~)2  -  K2Tv2] 

*  b*  1 


(eq.  3-33) 


As  before,  P^Fp  is  expressed  as  a  function  of  FQ/t^. 

The  pressure,  P,  is  estimated  from  the  dynamic  pressure, 

*  -  »vuv2' 


P  o 
Mv  v 


P  u 
Ks  V 


[F0/t“  -  /ir  In  Tg/Tv] 


nFp/t;)2  -  LIvj 


(eq.  3-34) 


In  the  case  of  long  pulses,  the  expression  for  pv/Fq 
is  different.  As  before. 


F  /t~  - 
O  P _ 

[L  +  C_T, 


s 


s  V' 


Here,  one  can  use 


(eq.  3-35) 


U. 


AZ 

<*P  '  V 


AZ 

'S 


</t~- 

P  ^ 


(eq.  3-36) 


Again  using  conservation  of  mass,  p  U.  =  p  _U  one  can 

V  V  s  s 

substitute  and  arrive  at  the  equation  for  pv/FC): 


p 


v 


F 

o 


K2T  2 

V^sV'V^  * 


-] 


o  p 


(eq.  3-37) 


The  pressure  is  found  in  the  same  manner  as  for  the  short 
pulse  and  is: 


_P 

F 


s-  »  U. 


_iVVL!s2iViL_ 

KJT  2 

IL  +  W  [Fe/t  -  — ^;1 


(eq.  3-38) 


7.  Limit  on  T_  to  Avoid  Ionization 
11  ■  . s - ...  — 

The  question  now  is  what  value  can  Tg  take  in  the 
elevated  surface  temperature  model.  As  pointed  out  previously 
if  the  temperature  is  allowed  to  become  on  the  order  of  2-3  eV 


i 
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then  the  ionization  degree  reaches  a  critical  value  where 
the  vapor  plume  begins  absorbing  the  laser's  radiation.  This 
absorption  by  the  plume  occurs  at  plasma  densities  given  by 
the  plasma  frequency 

2 

a)  =  n  e  /e  in  (eq.  3-39) 

p  e  o  o 

where  mQ  is  the  free  electron  density  in  the  plume.  If 
electron  densities  are  allowed  to  rise  above  this  cut-off 
value  for  propagation  of  electromagnetic  waves,  then  strong 
absorption  will  occur  in  the  low-temperature  plasma  plume. 

This  condition  for  cut-off  can  be  rewritten  as 

ne  >  1.12  x  1013/Ap2  (eq.  3-40) 

where  Ap  is  now  the  laser  wavelength  in  centimeters  and  ng 
is  the  electron  density  in  electrons  per  cubic  centimeters. 

If  we  assume  one  free  electron  per  atom  and  solid  density  in 
the  vapor  at  the  target  surface  then  Saha's  equation,  [23], 
can  be  used  to  predict  the  maximum  allowable  surface  temper¬ 
ature  to  avoid  cut-off  degrees  of  ionization 

2®  .  2.4  x  1015  I  .-IWOOOi/T,  (eq.  3_41> 

"n  8 

where 

n  =  neutral  particle  density 


U.  =  ionization  potential  of  the  target  material 
in  eV. 

The  degree  of  ionization  is  given  by  n  /n.  where  n,  *  n  +  n 

■*  e  t  ten 

and  can  be  written  as  follows: 


2.4  x  lo 


-15 


3/2 


-11600U./T 
i  s 


and 


n 


n. 


-C  + 


/c2  +  4n.  C 


2n, 


(eq.  3 


We  now  have  a  rough  temperature  criterion  for  significant 
plasma  production  to  occur.  Namely,  arrange  for  the  value 
PQ/t~  such  that  Tg  will  always  be  below  the  value  required 
for  cut-off. 


-42) 


n 


<  1.12  x  1013/X  2 


IV.  RADIAL  EXPANSION  OF  THE  VAPOR  PLUME 


As  the  vapor  plume  propagates  along  the  axis  it  will 
also  expand  radially.  To  assess  the  feasibility  of  the 
approach  one  needs  to  know  whether  the  radial  expansion  is 
commensurate  with  the  expected  geometry.  Therefore  we  need 
an  estimate  of  the  radial  expansion  as  a  function  of  the 
distance  from  the  target  surface. 

A  simple  estimate  can  be  obtained  if  we  assume  that  the 
speed  with  which  the  plume  radius  expands  is  equal  to  the 
sound  speed,  see  Figure  16. 

We  consider  therefore  a  slice  of  the  plume  which  has 
traveled  along  the  axis  for  a  time,  t.  For  the  purpose  of 
this  discussion  we  neglect  any  expansion  of  the  slice  in  the 
axial  direction  due  to  possible  pressure  gradients  in  that 
direction.  As  the  slice  propagates  along  the  axis  and 
expands  radially  the  temperature  and  density  must  change  and 
so  the  sound  speed  is  a  function  of  the  temperature  which  is 
a  function  of  the  z-position  on  the  axis. 


dR(t) 

”3t 


a(T) 


/cp(Y-l)T 


(eq.  4-1) 


The  expansion  is  assumed  to  be  adiabatic  so  that 


(eq.  4-2) 
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dR/dt  ■  a 


1 


dR/dt 


FIGURE  16.  MODEL  ILLUSTRATING  THE  RADIAL  EXPANSION 
OF  THE  VAPOR 


where  T  ,  P^  and  are  measured  at  the  target  surface  and 
o  o  o 

T^,  P^,  and  p^,  are  taken  to  be  at  some  later  time,  after 
the  vapor  has  expanded.  The  quantity  y  is  the  ratio  of 
specific  heats  and  is  1.66  for  monatomic  metal  vapor. 

For  small  changes  in  density  and  temperature,  dp  and  dt. 


(P  +  dpjV 

p 


,T  +  dT»  y/ y-1 

l  ip  * 


(eq.  4-3) 


+  dp  /T+_dTv  1/y-l 

p  '  T  ' 


(eq.  4-4) 


Using  binomial  expansion  for  the  temperature  term  yields 


l  +  2fi. 

P 


t  .  ,  1  v  dT 

1  +  (Tn:)  t" 


(eq.  4-5) 


or  simply. 


(1,  « 
'7-T  T 


(eq.  4-6) 


Since  the  density  decreases  with  increasing  radius  (again 
assuming  that  the  thickness  of  each  slice  does  not  change) 
one  can  write 


(p  +  dp)  (R  +  dR) 


p  R2  , 


(eq.  4-7) 


and 


(1  +&-)  (1  +  ^)2 
P  R 


(eq.  4-8) 


Again  using  binomial  expansion  for  the  quadratic  in  R  and 
substituting  for  dp/p  the  above  expression  yields 


(1  +  X>  =  1 


(eq.  4-9) 


By  assuming  the  second  order  (dT) (dR)  term  to  be  negligible 
compared  to  the  other  terms  and  rewritting  gives. 


1  dT 
y— 1  T 


(eq.  4-10) 


Rq  is  now  taken  to  be  the  initial  plume  radius  which  is 

the  laser's  focal  radius.  T  is  assumed  to  be  the  surface 

o 

temperature.  Integrating  both  sides  yields 


v 


o 


(eq.  4-11) 


(i)1"* 

V 


(f)1/2 

o 


(eq.  4-12) 


Substituting  /T  ■  a//Cp  (y — i )  and  a  *  dR/dt  yields 


( JL) 1-Y 

V 


dR/dt 


(eq.  4-13) 
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The  term  /Cp (y-1) Tq  is  simply  Uv,  the  initial  vapor  velocity. 
Therefore , 


and/ 


.  R  .  1— y  1  dR 

V  *  v;  « 


dR 

?=7 


Integrating  both  sides: 


t 

/  dt 

t 


or. 


1  + 


) 


(eg.  4-14) 


(eg.  4-15) 


(eg.  4-16) 


(eg.  4-17) 


But,  Uv  *  (t  -  t  )  is  the  axial  distance,  i,  that  the  vapor 


plume  will  travel  in  the  cavity.  Therefore, 


This  expression  can  now  be  used  to  calculate  the  radial 
expansion  as  a  function  of  the  axial  distance  traveled  i. 
The  previously  derived  equation  for  temperature  can  be 
used  to  calculate  T  at  that  distance  R: 


(JL)-2 

V 


The  density  and  pressure 
expansion  relationship: 


R  , 

(— )2 
'  R  ’ 


(Tji/Y-i  (eq.  4-20) 

Ao 

can  be  obtained  using  the  adiabatic 


l/Y-l 


(eq.  4-21) 


It  must  be  noted  again  that  this  calculation  gives  merely 
an  estimate  for  I  and  in  a  more  rigorous  model  would  have 
to  be  replaced  by  a  proper  fluid  dynamical  model. 


i 


0 
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V.  RESULTS 


There  are  now  two  cases  of  the  mass  removal  process. 
Both  involve  Fq  low  enough  to  avoid  significant  ionization 
so  that  all  the  laser's  energy  is  either  absorbed  at  the 
material's  surface  or  is  reflected.  There  will  be  little 
absorption  in  the  vapor  plume. 


A.  CASE  I.  SMALL  PULSE  TIME 

This  is  the  situation  where  the  following  condition  is 
is  met:  2  2 


Fo  > 


F  /t“  -  p  C  T  /k 

F  /t~r  °  P  — 3  3  -v _ 

op1  Pe<  L  +  CJO 


s  v* 


In  this  case  the  surface  temperature  will  be  elevated  slightly 
above  the  vaporization  temperature  and  is  described  by  equation 
(3-28) .  Figure  17  is  a  plot  of  the  elevated  surface  tempera¬ 
ture,  Ts  as  a  function  of  FQ/t^,  the  controlling  parameter 

at  10  urn.  On  this  plot,  the  lowest  value  of  T  for  each 

s 

material  corresponds  to  Tv  and  FQ/t^.  The  maximum  value 
for  Tg  is  the  temperature  at  which  there  will  be  cut-off  and 
significant  absorption  of  the  radiation  occurs  by  the  vapor 
plume  due  to  ionization.  That  maximum  Tg  was  calculated  by 
first  finding  the  maximum  allowed  free  electron  density,  ne, 
given  by  equation  (3-40) ,  Next,  calculate  the  maximum 
allowed  degree  of  ionization,  ne/nf  Then,  use  Saha's 


equation,  equation  (2-42)  to  find  the  maximum  allowed  T  . 

9 

It  will  be  a  different  value  for  each  material  and  wave¬ 
length,  and  can  vary  by  several  thousands  of  degrees. 

Figure  18  is  a  plot  of  MFQ  vs.  FQ/t^  for  Case  I  of 
short  laser  pulses  and  Case  II  of  long  time  pulses  from 
equations  3-20  and  3-21.  The  difference  between  Case  I  and 
Case  II  is  slight  except  at  high  values  of  FQ/t^.  Case  I 
gives  a  lower  value  of  MFq  than  Case  II.  In  this  plot,  the 

maximum  value  for  Case  I  corresponds  to  a  F  /t~  and  a  T 

op  s 

where  plume  absorption  begins.  Case  II  was  also  plotted 
to  that  same  value  of  F  /t~,  however,  the  surface  tempera- 
ture  does  not  rise  above  vaporization  temperature  and  could 
be  plotted  to  much  higher  values  of  FQ/t^. 

B.  CASE  II.  LONG  PULSE  TIMES 

Case  II  is  the  situation  where  the  condition  on 


Fo  > 


v2_  2 

U  [  (F  /t")2 - 2-] 

V  O  P  _ K 

F  /tT  -  p  C  T  Sk 

r0/V  °  »  s  * v 


p.lL  +  c.V 


is  not  met.  This  is  the  case  of  low  F  and  long  t  .  The 

o  P 

surface  temperature  remains  at  essentially  Tv>  In  this 
case  abosrption  by  the  plume  never  occurs . 


C.  MASS  REMOVAL  CALCULATION 

A  numerical  example  is  now  presented  to  show  the  use  of 
the  various  equations  and  to  illustrate  that  masses  on  the 


53 


Ag(II) 


_  ^4  ^  -2  ®  1/2  ® 

/t  x  10  Watts -an  -sec 
5  P 


F  /t  FOR  CASE  I,  SHORT  PULSE 
O  p 


[•] 


order  of  milligrams  could,  theoretically,  be  removed  without 
plasma  production  at  a  materials  surface.  The  material 
chosen  for  the  example  is  Zn,  as  it  has  a  relatively  high 
Z-number  and  low  latent  heats  and  vaporization  temperature. 
The  following  thermal  properties  are  from  Table  II: 

p  =  7.13  g/cm^ 

s 

C  =0.39  Joules/g°K 
s 

2 

<  «  0.41  cm  /sec 

K  *  1.13  Watts/cm  °K 

Tv  »  907  °C 

Z  *  65.38 

Let  us  assume  that  F^/tT  *  10*  Watts  sec  '  /cm  and 

o  p 

-1  2 

that  the  laser's  focal  area  *  10  cm  .  First  the  condition 
for  the  two  cases  must  be  tested  from  equation  3-28,  Figure 
19. 


F  *  6.22  x  io8  Watts/cm2 

o 
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g 

If  F  is  chosen  to  be  less  than  6.22  *  10  then  Case  II 
o 

must  be  used.  If  not,  Case  I  must  be  used.  Let  us  choose 
Fq  -  2  x  io®  Watts/cm2.  Then, 

t  -  (loVFi2  =  2.5  x  10~5 

P  ° 

=  25  ysec. 

From  equation  3-21  or  Figure  18  for  Case  II, 

4  -4 

MF  *  3.79  x  10  Watts  grams  cm 

o 


or. 


M 


3.79  x  10 


-2 


grams  cm 


-2 


-1  2 

If  the  focal  area  is  assumed  to  be  10  cm  ,  then  the  mass 
removed  is  3.79  mg.  The  axial  expansion  during  the  pulse 


l  •  U(t-t)  *  0.83  cm. 

s  p  V 

The  density  and  pressure  before  radial  expansion  can  be  found 
using  equations  3-37  and  3-38. 


Pv  -  2.29  x  10“2  g/cm3 


7  2 

2.65  x  10  ergs/ cm 


3.85  psi 


F  /t  x  104  Watts-cnf  2-sec1/ 
o 

FIGURE  19.  CONDITION  FOR  CASE  I  OR  CASE  II,  EQUATION  2-28 


If  the  anode-cathode  separation  is  assumed  to  be  4  cm  and 
Rq  »  0.18  cm,  then  the  radial  expansion  at  z  ■  4  cm  can 
be  found  from  equation  4-19. 

R  *  1.65  cm. 

The  pressure,  density  and  temperature  after  radial  expansion 
to  R  *  1.65  cm  and  axial  expansion  to  z  *  4  cm  are  calculated 
from  equation  4-21. 


P  =  2.39  x  10“3  psi 

p  *  2.73  x  10“4  g/cm3 

T  *  61.5  °K 

This  shows  that  a  relatively  low  power  C02  laser  pulse 
6  2 

of  2  x  io  Watts/cm  of  25  psec  duration  can  produce  a 
zinc  vapor  plume  with  the  desired  mass  and  dimensions. 
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VI.  CONCLUSION 


As  stated  in  the  Introduction,  the  purpose  of  this  study 
was  to  assess  the  general  feasibility  of  using  a  laser's 
radiation  to  remove  masses  on  the  order  of  milligrams  from 
a  target's  surface.  This  material  will  fill  a  discharge 
channel  between  the  anode  and  cathode  of  a  high  voltage 
pulse  device  which,  when  activated,  will  couple  electrical 
energy  into  the  material  in  the  discharge  channel  causing 
the  atoms  to  undergo  electronic  transitions  and  then  giving 
off  intense  x-rays. 

The  model  presented  herein  is  a  heat  balance  model  at 

the  target  surface.  No  attempt  has  been  made  to  solve  the 

complicated  heat  conduction  equations  with  its  associated 

moving  boundaries.  In  this  model,  the  controlling  parameter 

is  FQ/t ,  the  absorbed  power  density  and  the  laser  pulse  time. 

Two  separate  cases  emerge  from  the  model.  Case  I  is  the  case 

where,  keeping  FQ/t  constant,  FQ  takes  on  a  relatively  high 

value  while  the  pulse  time  is  small.  Case  II  considers 

the  situation  of  the  same  value  for  F„/t  but  smaller  F_ 

o  o 

and  longer  pulse  time.  Using  this  model,  theoretically, 
masses  of  the  order  of  milligrams  can  be  removed  from 
metallic  surfaces  with  relatively  low  power  densities. 

The  mass  removal  model  also  predicts  low  temperatures 
at  the  target  surface  and  thus  in  the  discharge  channel. 

As  stated  in  section  II  this  will  produce  low  degrees  of 


ionization.  In  order  to  increase  the  coupling  efficiency 
of  the  electrical  energy  into  the  vapor  plume ,  the  ioniza¬ 
tion  level  must  be  increased.  This  could  be  done  by  increasing 
the  laser  power  level  towards  the  end  of  the  pulse  or  by  using 
microwave  radiation  after  the  plume  has  propagated  into  the 
discharge  channel  and  prior  to  the  application  of  the  high 
voltage  pulse. 

In  addition  to  mass  removal,  the  controlling  parameter 
F  /tT  is  used  to  predict  the  temperature,  density  and 
pressure  of  the  vapor  plume  as  it  begins  to  propagate. 

As  the  plume  propagates  into  the  vacuum  chamber  it  undergoes 
axial  and  radial  expansion.  Temperature,  density  and  pres¬ 
sure  decrease  during  this  process  and  can  be  predicted  using 
a  model  of  adiabatic  jet  expansion  presented  in  section  II. 

The  mass  removal  model  predicts  that  a  laser  can  remove 
masses  from  a  target  surface  large  enough  to  be  used  in 
the  proposed  experiment.  In  addition,  the  expansion  of  the 
vapor  plume  is  sufficiently  small  as  not  to  damage  the  com¬ 
ponents  located  in  the  vacuum  chamber.  Using  a  long  laser 
pulse  (Case  II)  the  proposed  technique  shows  that  the  vapor 
in  the  discharge  channel  is  relatively  cold  and  requires 
pre-ionization  prior  to  application  of  the  high  voltage 
pulse. 
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APPENDIX  A 


The  interactions  of  materials  with  electromagnetic 
waves  is  formulated  on  the  basis  of  Maxwell's  equation 


7  x  E 


3H 


3M 


U0  3t  W0  3 t 


—  —  3E  3P  — 

7  x  H  "  eo  if  +  af  +  J 


7  «E  =  -  -i-  7  •  P  +  p 

eo  e 

7-h  *  -  7*m 

where, 

E  =  Electric  field 

H  s  Magnetic  field 

M  =  Magnetization 

J  =  Current  density 

P  =  Polarization 

e  =  Dielectric  constant 
o 

uQ  =  Magnetic  permeability 

p  =  Charge  density 
© 
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We  will  assume  that  our  concern  is  with  non-magnetic , 
electrically  neutral  materials.  In  conducting  mediums  it 
is  the  contribution  from  the  free  electrons  that  predominate. 
This  will  allow  us  to  ignore  the  polarization  effects. 
Maxwell's  equations  now  reduce  to  the  following: 


7  x  E 


U 


3H 
o  dt 


7  X  H  *  EoS+J 
7-E  -  0 

7-H  -  0 


The  general  wave  equation  for  E  is  then  found  by  taking  the 
curl  of  the  first  equation  and  the  time  derivative  of  the 
second  equation. 

V  x  V  x  E  -  ?  x<-  „0  ||) 


7ilf  x  ff) 


l_(e  »)  +  M 

3tkeo  3t'  3t 


Using  the  vector  identity 


7  x  7  x  E  ■  7(7*E)  -  72  E 
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and  noting  that  v.E  =  0,  we  can  exchange  the  order  of 
differentiation  in  the  electric  field  equation  and  multiply 
the  magnetic  field  equation  by  yQ  to  give 


?2  5  "  "o  lt(7  *  H) 


3/77  TTi  .  32E  .  3  J 

-  »o  at'7  -  H»  +  Vo 


2  1 

Adding  and  substituting  C  * 


woeo 


_2  F  1  3E  3 J 

7  E  ’  ^  3t  =  wo  at 


Plane  harmonic  waves  solutions  for  the  electric  field  are 
of  the  form 

E  -  Eq  exp{  [-2irkz/A  J  +  tiu  (t  -  nz/C)  ] } 

where 

n  =  Real  part  of  index  of  refraction 
k  s  imaginary  part  of  index  of  refraction. 


Substituting  this  into  the  electric  field  equation  derived 
above  along  with  Ohm's  law,  J  ■  oE  -  -Nev,  gives  us  an 
expression  for  n-ik. 


o  is  the  conductance  and  is  in  general  a  complex  quantity. 


namely  o  »  o^-ioj*  Then, 


n  -  ik 


e  w 

o  o 


In  order  to  evaluate  n  and  k  one  can  use  the  Drude-Lorentz 
model.  In  this  model  the  atom  and  its  outermost  electron 
are  assumed  to  follow  the  equation  governing  the  motion  of 
the  harmonic  oscillator: 


m* 


dv  m*v 
dt  t 


-  e  E 


where 


m*  =  effective  mass 
t  =  relaxation  time. 

If  a  ^elution  of  the  form  e_1“t  is  assumed  for  v  and  E, 
and  Ohm's  law,  oE  *  -Nev  is  substituted  for  v,  one  can  show 
that 


Ne  t 

0  “  m*  ( 1  +  u  t  ) 

By  separation  of  the  real  and  imaginary  parts  one  can 
write  an  expression  for  and  02* 


o  U)T 

O 

"  T7 

1  +  U  T 


where  aQ  is  the  static  of  dc  conductivity.  As  the  frequency 
approaches  zero,  u  -*■  0, 

Ne2T 

°o  *  m* 


The  real  and  imaginary  parts  of  the  index  of  refraction, 
n  and  k  can  now  be  found  using  a  complex  identity: 


2n 


2 


+ 


+ 


x  can  now  be  eliminated  from  the  expression  for  o^,  and 
by  using  the  equation  for  aQ.  The  resulting  equations  for 
and  a 2  are: 


1  +  u2  ( - £)  2 


Ne 


2  ro*°o 


Ne 


0  ra*a«  2 
1  +  u2  ( - %)2 


Ne 
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If  normal  incidence  is  assumed,  the  abosrbtivity  of  a 
material  can  be  calculated,  again  using  Maxwell's  equations 
and  the  appropriate  boundary  conditions.  Figure  20  illus¬ 
trates  the  geometry. 


V  x  E 


-M 


3 1 


V  x  H 


a  E  + 
x 


As  before. 


E  *  Eq  exp{  [-2tt1cz/A  ]  +  [iiu  (t  -  nz/c)  ]  } 

Substituting  this  expression  for  electric  field  into 
Maxwell's  equations  yields  a  relationship  between  H  and 
E: 


MC 


The  boundary  conditions  at  the  material  interface  are: 


The  subscripts  1  and  2  are  for  the  mediums.  Medium  1  is 
vacuum  and  medium  2  is  the  target  material.  The  subscripts 
i,  r,  and  t  denote  incident,  reflective  and  transmitted 
waves.  Using  these  boundary  conditions  and  the  relationship 
for  Hy  in  terms  of  Ex,  the  ratio  of  E^/E^.  can  be  readily 
obtained: 


n  -  1  -  ik 
n  -  1  +  ik 


r 


where  r  is  the  amplitude  reflection  coefficient.  It  has 
been  assumed  that  n^  *  1  and  *  0,  which  is  the  case  in 
vacuum.  The  reflectivity,  R  is  given  by  |r  *  r|  and  there 
fore  is  written  as: 


n  ,  (n-1) 2  +  k2 

(n+1) 2  +  k2 


The  absorbtivity.  A,  is  the  desired  result  of  the  derivation 
where : 


A  =  (1  -  R) 

The  abosrbtivity  can  be  written  as: 

A  = 

(n+1) z  +  k 
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